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ABSTRACT 

We consider the effect of free-streaming axion emission on numerical 
models for the cooling of the newly born neutron star associated with 
SN1987A. We find that for an axion mass of greater than ~10 -3 eV, axion 
emission shortens the duration of the expected neutrino burst so significantly 
that it would be inconsistent with the neutrino observations made by the 
Kamiokande II (KII) and Irvine-Michigan-Brookhaven (IMB) detectors. 
However, we have not investigated the possibility that axion trapping (which 
should occur for masses ^0.02 eV) sufficiently reduces axion emission so that 
axion masses greater than ~2 eV would be consistent with the neutrino 
observations. 
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I. INTRODUCTION 

SN1987A confirmed astrophysicists' most cherished beliefs about the nature of type 
II supernovae: that they are associated with the formation of neutron stars and that in the 
process they release their binding energy in thermal neutrinos. 1 In addition, it has also 
provided a wealth of information about the properties of neutrinos and other hypothetical, 
weakly-interacting particles. In particular, the detection of 11 neutrino events over M2 sec 

by the Kamiokande II (KII) detector^ and 8 neutrino events over ~6 sec by the Irvine- 
Michigan-Brookhaven (IMB) detector 3 indicates that thermal neutrinos with temperature -4 
MeV indeed carried away the bulk of the ~(2-4)xI0 53 ergs of binding energy from the 
explosion. 1 In turn, these observations have led to constraints to the mass, charge, unknown 
interactions, magnetic moment, speed of propagation, and lifetime of the electron antineutrino, 
on the possible existence of right-handed neutrinos and their coupling strength, and on the 

possible existence and mass of the axiom 4 -* It is the last of these issues which we will 
address in this paper. 

The axion was proposed in 1977 to solve the strong CP problem of quantum 
chromodynamics (QCD). 9 To date, it is still the most attractive solution to this solitary 

blemish on QCD. The axion necessarily couples to nucleons, with a strength proportional to 
its mass, and may also couple to electrons (though it need not). Astrophysical arguments 
(red giant emission) preclude an axion of mass greater than -0.01 eV (DFS-type axion) 10 or 
-3-60 eV (hadronic-type axion), 11 while the cosmological production of axions precludes an 
axion of mass less than a few x 10* 6 eV.'2 Thus, there exists a window of allowed axion 
masses: few x 10 -6 eV - (3-60) eV (hadronic) or few x 10 -6 eV - 10 -2 eY (DFS). 

For axion masses in this window, axion emission by nucleon-nucleon axion 
bremsstrahlung from the newly born neutron star associated with SNI987A should have been 
a significant cooling mechanism. If the axion exists, such a heat sink would have 
accelerated the cooling and thereby shortened the duration of the neutrino signal that was 
detected by the KII and IMB detectors. With neutrino cooling alone, theoretical protoneutron 
star models predicted that the neutrino burst would last for on the order of several to many 
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seconds 13 — not inconsistent with the observations. Several groups of authors 4-8 have 
argued that an axion with mass in the range (2.0 - 10~ 3 ) eV is ruled out, as, for such a 
mass, axion emission would be so important that it would drastically reduce the duration of 
the neutrino burst. [For axion masses less than ^0.02 eV, axions freely stream out of the 
core; for masses greater than ~0.02 eV, axions become trapped and are radiated from an 
M axion sphere;' 14 and for axion masses greater than -2 eV, trapping apparently renders axion 
emission insignificant. 4 ] 

The purpose of the present work is to address in a very quantitative way the 
axion mass limit in the free-streaming regime (m a < 0.02 eV); in a later work, we plan to 
address the trapped regime (m a > 0.02 eV). In particular, while previous authors have used 
axion emission rates that are valid in either the strongly degenerate regime or the 
nondegenerate regime, here we use the exact rates computed by Brinkmann and Turner. 5 
Previous authors have either neglected the back reaction of axion cooling on the model of 
the cooling neutron star or tried to incorporate axion cooling in an ad hoc manner. In this 
work we fully and self-consistently incorporate axion emission into the model for the cooling 
of the young neutron star. To derive a limit to the axion mass, we use the expected 
duration of the delected neutrino burst in the KI1 and 1MB detectors, whereas in the 
previous works, the axion luminosity or the total energy emitted in axions, neither of which 
were observable, were used. Because of uncertainties in both the equation of state at 
supranuclear densities and the actual baryon mass of the remnant, we explore a variety of 
numerical models 14 to test the sensitivity of our limit to the theoretical models of the nascent 
neutron star employed. 

As mentioned above and as expected, the observable most sensitive to axion 
cooling is the duration of the neutrino burst. For the wide range of models we have 
explored, axion emission has virtually no effect on the burst duration if the axion mass is 
less than or equal to 10' 4 eV. However, if the axion mass is greater than or equal to 10~ 2 
eV, the duration of the neutrino bursts in both detectors for all models considered is less 
than 1 second. Such a short time is clearly in conflict with the observations. For all the 
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models we considered, the neutrino burst duration dropped percipitously at an axion mass of 
10 -3 eV. strongly suggesting that the upper limit to the axion mass (in the free-streaming 
regime) is 10 -3 eV (to within approximately a factor of 2). 

The paper is organized as follows: in Sec, II we describe the numerical models 

we use to simulate neutron star cooling and to compute the expected neutrino signals in the 
KII and IMB detectors; in Sec. Ill we discuss the axion and the axion emission rates we 
use; in Sec. IV we describe the effect of axion emission on the cooling of the nascent 
neutron star associated with SN1987A and on the expected neutrino signals; we end with 
discussion and a summary in Sec. V. 

II. DESCRIPTION OF THE NUMERICAL CODE AND MODELS 

For the purposes of this study, the protoneutron star evolution code of Burrows 
and Lattimer 13 and Burrows 14 was modified to include the axion energy loss rates recently 
derived by Brinkmann and Turner 5 and described below in Section III. The code uses 
standard relaxation techniques to solve the general relativistic equations of stellar structure. 
It incorporates all relevant redshift factors, follows all six neutrino species (three two- 
component neutrino and antineutrino species), employs a "realistic" nuclear equation of state 
(EOS), 13 - 14 and has a sophisticated neutrino opacity algorithm. The neutrinos are assumed to 
be thermalized with the local matter temperature and to be emitted with a Fermi-Dirac 
energy distribution. 

The core of a massive star becomes unstable upon reaching the Chandrasekhar 
ma<;<; (~1.4 M 0 ) and implodes. Core collapse proceeds through five orders of magnitude in 
central density and two orders of magnitude in radius, and is halted only when the matter 
stiffens upon reaching nuclear densities. The inner core rebounds into the outer core, a 
shock wave is formed, and the inner structure, the protoneutron star, rapidly achieves 
hydrostatic equilibrium. It is during the quasi-hydrostatic neutronization and cooling phase 
(timescale ~ seconds) of the protoneutron star, not during the dynamical phase of collapse 
and shock wave formation (timescale ~ milliseconds) that the prodigious neutron star 
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gravitational binding energy (~2-4x 1 0 53 ergs) is released. In the standard model, the energy 
is radiated as neutrinos (and antineutrinos) of all species. These protoneutron star neutrinos 
constitute the signature of "core collapse" (cf.. Ref. 14) that was apparently detected by the 
KII 2 and 1MB 3 detectors. At the high densities and temperatures typical of a nascent 

neutron star, even neutrino mean-free-paths (\„) are small compared to the size of the core. 
Therefore, neutrino cooling proceeds on a long "diffusion" timescale (seconds) and not on the 
short production or light-travel timescales («second). Hence, the neutrino signal is spread 
over many seconds. The neutrino signal can be separated into two phases. The first is an 
outer mantle cooling phase, powered in part by residual accretion and quasi-static mantle 
collapse during the first 1-2 seconds. The second phase is a later, longer (> 2 seconds) 
phase of inner core cooling during which the neutrino luminosity is powered by neutrino 
transport of energy from the core, characterized by the longer neutrino diffusion timescale 
(several sec). During both phases, the neutrinos escape from the periphery at the 

"neutrinosphere" where ~ R (radius of the neutron star ~ 10 km). Because the capture 

process (t> e + p -» n + e + ) has a much larger cross section than the various scattering 

processes ( v l + e - ■+ v % + e”), the neutrino signal is dominated in H 2 0 detectors by its H e 
component. In H 2 0 Cherenkov detectors, the secondary positron (or electron for a scattering 
process) is detected by its Cherenkov light. As has been pointed out by many analyses 
now, 1 the IMB and KII detections are consistent with an effective v e temperature (T. ) of 
~4.0 MeV, a characteristic cooling time of ~4 seconds, and a total binding energy (6 x E. ) 
of ~2-3xl0 53 ergs, all consistent with the standard model of a Type II supernova. However, 
if the axion exists and has a large mass (m a ), and therefore a large coupling, the resulting 
axion energy losses would be at the expense of neutrino emission, and the signal in neutrino 
detectors would be altered. In this paper, we investigate the implications of axion emission 
for the predicted signals in IMB and KII. Heretofore, authors have used the results of 
previous numerical models (without axion emission ), i.e., temperature and density profiles, to 
compute axion emission. However, to properly take account of the feedback on the 
temperature-dependent axion emission of the axion-cooling induced temperature decreases, a 
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detailed evolutionary calculation, which incorporates axion emission ab initio is required. By 
doing so. we are able to sensibly address the question: For what range of axion masses 

could we not fit the SN1987A neutrino data? 

Three protoneutron star models from the more comprehensive work of Burrows 14 
were evolved at six different axion masses between O.OeV and 10~ 2 eV for 20 "physical" 
seconds after bounce. In addition, it was assumed that axions once produced freely stream 
out — a valid assumption for m a < 0.02 eV. 4 Using the published detector fiducial masses, 
efficiencies and energy thresholds, the appropriate neutrino interaction cross sections, and a 
supernova distance of 50 kpc, the predicted neutrino signals in both KII and IMB from 
SN1987A were calculated. The effect on this signal of axion emission will be described in 
Section IV. The models (A, B, and C) were chosen to represent a range of possibilities, 
since the precise neutron star EOS and the ultimate baryon mass (Mg) of the residue are not 
accurately known. The initial entropy and lepton profiles employed were similar to those 
found in the collapse literature. 15 Model A is model 57 from Ref. 14. which starts at Mn - 
1.3 M q and accretes to a large 1.8 M 0 by means of an accretion rate that is taken to decay 
exponentially with a time constant of 0.5 seconds. The EOS employed in Model A is stiff, 
as described in Ref. 14. Similarly, Model B is stiff model 55 of Ref. 14, in which an 
initial core of mass 1.3 Mq accretes to a mass of 1.5 Mq with a similar accretion time 
constant. Model C is soft model 62 from Ref. 14 that in all ways, save stiffness, is the 
same as Model B. A soft calculation with baryon mass parameters similar to those of model 
A is not included in the set because the maximum baryon mass for the soft EOS is so small 
(1.6 M 0 ) that a black hole would form early on (<1 second), truncating the neutrino 
emission. With Models A, B, and C, we represent a range of realistic behavior and binding 
energies of the protoneutron star. 

III. THE AXION AND AXION EMISSION RATES 

The axion is the hypothetical pseudo Nambu-Goldstone boson associated with the 




where N is the color anomaly of the PQ symmetry. 16 Generically, there are two types of 
axions: axions that couple to both quarks and leptons, with strength ^m/f a (m - quark or 

lepton mass), the so-called DFS-type axion; 17 and axions that only couple to quarks, and 
perhaps not even to the ordinary light quarks, but only to heavy, exotic quarks, the so-called 
hadronic-type axion. 18 Both types of axions couple to photons and nucleons through 
electromagnetic and color anomalies. 

The relevant couplings of both types of axions to electrons, photons, and nucleons 
are summarized in Ref. 4, and discussed in detail in Ref. 16. For the purposes at hand, 
we are only interested in the axion-nucleon couplings, as by far the dominant axion emission 
process from SN1987A is nucleon-nucleon, axion bremsstrahlung. 8 Those couplings, as 
computed in the naive quark model, are 4 


g an - [(X^/N - X^/4N) - 0.20] m/(f a /N) 

- 1.5 x 10- 7 [(XJj/N - Xy/4N) - 0.20] (m a /eV) 


(2a) 


gap « [(Xy/N - X|,/4N) - 0.55] m/(f a /N) 

- 1.5 x 10- 7 [(X^/N - X^/4N) - 0.55] (m a /eV) 


(2b) 


where X^ and X' u are the PQ charges of the up and down quarks, m is the nucleon mass, 
and the interaction Lagrangian (with nucleons) is 
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tf- mt “ • • • + (g an /2m) (n 7 m 7 5 n) + (g ap /2m) (p t m 7 5 P) 3* a . (3) 


The axion-nucleon coupling arises in roughly equal amounts from two sources: 
the direct coupling of the axion to up and down quarks (reflected in X‘ u and X' d ) and 
axion-pion mixing. For this reason, the axion-nucleon coupling is of the order of m/(f a /N), 
whether the axion is of the hadronic or of the DFS type. For comparison, the hadronic 
ax ion-electron coupling, which arises due only to radiative corrections, is some four orders 
of magnitude smaller than that of a DFS axion. 

The Feynman diagrams for nucleon-nucleon axion bremsstrahlung are shown in 
Fig. 1. The matrix element squared for this process has been computed by Brinkmann and 
Turner, 5 by H.-S. Kang, 19 and in the degenerate limit for the process nn - nn + a by 
Iwamoto. 20 In the nonrelativistic limit (i.e., to lowest order in T/m ~ few x 0.01), the 

matrix element squared is constant and is given by 5 


L i-«i 2 - ¥ s ^r L 


(3 - p) 


spin 


for the process nn - nn + a, by 


\MV- 


256 gl^f'nr 


m: 


(3 - fi) 


spin 


for the process pp pp + a, and by 

22 \m\ 2 - 256 f ‘ m! 


3mJ 


San San 

2 


(4) 


(5) 


(2 - 4/3/3) + 25 3~ - 4 — 


gan ~ t ~ gap 


(5 - 2/3/3) (6) 


spin 


for the process np - np + a. Here f - 1.05 is the neutral pion-neutron dimensionless 
coupling, and p = 3 <(2 f) 2 ) is the phase-space weighted average of the spatial dot product 
between the direction of the 3-momentum transfer in the direct and exchange diagrams, in 
the degenerate limit (p|/2m > 3T, or T ^ 20 MeV) 0-0, while in the nondegenerate limit. 
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/3 -* 1.0845. For complete details of the calculation of the matrix element squared, see 
Ref. 5. 

The axion emission rate (energy per volume per time) is given by a 15-dimensional 
phase-space integral: 


dn, dn 2 dlij dfl 4 dll a (2?r) 4 S 


x ^ \M I 2 5W(P! + p 2 - p 3 - p 4 - P a ) E a f, f 2 ( I -f 3 ) (l-f 4 ) , (7) 

spin 


where dll; - d 3 p i /(2rr) 3 2E; is the Lorentz-invariant phase-space element, the labels i - 1-4 
denote the incoming (1,2) and outgoing (3,4) nucleons, the label i - a denotes the axion, and 
S is the usual symmetry factor for identical particles in the initial and final states (S - 1 for 
np -» np + a; S = ‘/ 2 x l / 2 = '/« for nn -+ nn + a, or pp -* pp + a). The nucleon phase- 
space distribution functions f, are given by f, = [exp(E,/T - ju,/T) + l] -1 . Under the 
assumption that the matrix element squared is constant (which is accurate to about 10-20%), 
Brinkmann and Turner 5 have evaluated this 15-dimensional integral numerically. Summing 
over all three bremsstrahlung processes, they find 


c? a - 64(m 2 - 5 T 6 ' 5 /mJ) f« [(1-/3/3) g 2 an I(y„ yi ) + (1-/J/3) g 2 p I(y 2 .y 2 ) 


+ m-M 
9 



+ Sap 
2 


i(yi.y 2 ) + 


4(6-4<3) 

9 



( 8 ) 


where the first term accounts for nn -*• nn + a, the second term pp -♦ pp + a, and the third 
and fourth for np -* np + a. The quantities y t - 4 n /T and y 2 - n p / T. where /x n and |i p 
are the chemical potentials of the neutron and proton. The quantity l(y !,y 2 ) is a three- 
dimensional integral that must and has been evaluated numerically. A convenient analytical 
expression (accurate to better than 25%) and a "look-up table" (accurate to better than 5%) 
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are given in Ref. 5. 

To compute <S a one must specify g an , g ap , and 0. The axion-nucleon couplings g an 
and g ap are obviously model-dependent, and /? depends upon the degree of degeneracy. For 
definiteness, as well as simplicity, we will take 

gan * gap = - 7.6 x 10' 8 (m a /eV) (9) 

and )3 - V 2 . Given the other uncertainties inherent to this problem, and the fact that any 

- 1/2 

axion mass limits derived scale only as <5 a , these simplifications seem well-justified. In 

any case, a more specific treatment is always possible. 

IV. THE EFFECT OF AXION COOLING ON THE SN1987A NEUTRINO SIGNAL 

The results from the 18 model calculations of this study ([models A, B, and C] x 

[m a (eV) * 0.0, 10"\ 3x1 0“\ 10~\ 3x1 0‘ 3 , 10’ 2 ]) are summarized in Figs. 2-6 and Table 1. 

Figure 2 depicts the dependence of both the total energy lost to all species of neutrinos (E„) 

and the total axion energy loss (E a ) as a function of m a for models A, B, and C. [Note: 

"total energy" here denotes the energy carried off in the first 20 seconds; by 20 sec, the 

total energy has essentially converged.] Though the calculations were performed for only six 
values of m a , continuous curves based upon interpolation are presented. As Fig. 2 
demonstrates, for low values of m a , E y falls in the range of reasonable neutron star binding 
energies and only gradually decreases with increasing m a . For all models, E u and E a begin 
to respond to increasing m a near ~3xl 0~ 4 eV, but only very gradually. Even for m a * 10“ 2 
eV, E a is only 45%, 56%, and 69% of E u at m a « 0.0 eV for models A, B, and C, 
respectively. This sluggish dependence of E a on m a is a consequence of the feedback of 
axion cooling on the axion losses themselves. The axion emission rate varies as p 2 T 3 * 5 
(nondegenerate limit) or p [ ^ T 6 (degenerate limit); and because of this temperature/density 
dependence, axion emission is most significant deep in the core, which, as mentioned earlier, 
holds only about V 2 the heat released in the formation of the neutron star. Because of the 


stiff temperature dependence, axion emission is self-quenching, and, further, because the core 
only contains -7 2 the heat, the total axion losses amount to only - l /2 the total binding 

energy. As a result, a factor-of-10 increase in m a from 10“ 3 eV to 10“ 2 eV, which without 
feedback would imply a factor-of-100 increase in axion energy losses (oc m 2 ), actually results 
in increase factors of 1.9, 2.8, and 1.6 for models A, B, and C, respectively. The lethargy 
of (E a , E^) vs. m a is also reflected in Fig. 3, which shows the slow decrease of N K and Nj, 
the expected event total in K1I and IMB, respectively, for all three models. By m a - 10“ 2 
eV, the predicted event totals (Nj s) have decreased less than 50%. Based upon E„ or the 
number of neutrino events, one would be hard-pressed to rule out an axion as massive as 
10’ 2 eV! 

However, as Fig. 4 indicates, the neutrino signal durations in the 1MB and KII 

detectors are sensitively decreasing functions of m a beyond -3x1 0“ 3 eV. In Fig. 4, At(90%), 
the time it takes the accumulated number of events to reach 90% of the final total number 
of events, is plotted as a function of m a [90% of the final total is an arbitrary choice; 
similar behavior would follow for the choice of 60% or 70%]. By m a - 10‘ 3 eV, At(90%) 
has plummeted to values inconsistent with the long duration of the KII and IMB detections, 
and for m a * 10~ 2 eV, the pulse duration for both detectors of all models is less than 1 sec. 
The cause of this can be traced as follows. The early phase of neutrino emission results 
from the initial heat in the outer mantle and accretion, and has a short timescale (-1 second) 
because both neutrino diffusion in the low density, outer mantle and residual accretion are 
rapid (<1 second). The residual heat of the inner core is transported by neutrino diffusion 
to the outer core and the neutrinosphere on a timescale of several seconds, and thereby 
powers the late time neutrino flux. If the inner core heat source did not compensate for 
the quicker outer core losses, the temperature of the neutrinosphere (e.g., T_ for t>'s) and 
the neutrino luminosities would dive after -1 second and the neutrino flux would shut off. 
For m a *= 0.0 eV, the early, short phase smoothy merges into the later, long phase that 

accounts for -V 2 of the signal. 14 However, the major effect of axion emission is the cooling 
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of this high-density inner core crucial to powering the second phase. As m a approaches 
10 -3 eV. the core is rapidly depleted of heat and cannot supply the energy for the second 

phase. To illustrate this effect, the time evolution of the matter temperature profiles (T(M)) 

for model A for m a = 0.0 eV and m a = 10 -2 eV are depicted in Figs. 5 and 6, respectively. 
Each curve represents a snapshot in time. The lowest curve is the initial profile. As the 
shock-puffed outer core settles from R - 10 2 kilometers to R « 10 kilometers, the resulting 
compression raises the temperature of the outer core dramatically. Subsequent accretion 
further compresses the protoneutron star, and temperatures near ~50 MeV are achieved. 
Such high temperatures are of course not manifested directly in the neutrino signals, since 
the neutrinosphere is located on the periphery where T ~ 3-5 MeV. In both Figs. 5 and 6, 
the snapshots are every 100 milliseconds for the first 2.0 seconds and then every 2.0 seconds 
until the end (t = 20.0 seconds). The large temperature spike in Fig. 5 drives a neutrino 

flux into the center, thereby raising its temperature and storing heat for phase two. In Fig. 

6, it is plain to see that efficient axion cooling has refrigerated the inner core completely 

and depleted the heat reservoir for phase two. 

The behavior of At(90%) in Fig. 4 echoes the above-described phenomenon. For 
instance, at m a = 10~ 2 eV in model A, T , instead of being ~4.0 MeV at t - 1.0 second, is 
a tepid ~2.5 MeV. In model A. at m a = 10 -3 eV, At(90%) for IMB is only 40% of its m a « 
0.0 eV value and at m a = 10 -2 eV, it is only ~13% of that value. Since the m a - 0.0 eV 
models fit the IMB and KII detections. Fig. 4 strongly suggests an upper limit to m a of 10 -3 
eV based on signal duration alone. 

V. DISCUSSION AND SUMMARY 

We have fully incorporated axion emission (in the free-streaming limit) into 

numerical models of the initial cooling of the newly born neutron star associated with 
SN1987A. The dominant process for axion emission is nucleon-nucleon, axion 

bremsstrahlung, and our rates for this process are taken from Ref. 5, where the matrix 
element for this process has been calculated exactly and the phase-space integrals have been 
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evaluated numerically. Based upon the predicted neutrino flux from our models and the 
published detector response parameters for the KII and IMB detectors, we have calculated 
the expected characteristics of the neutrino pulses which would have been detected for axion 
masses of 0, 10" 4 eV, 3x1 0“ 4 eV, 10~ 3 eV, 3x1 O’ 3 eV, and 10~ 2 eV. By comparing the 
expected characteristics of the neutrino pulses with the neutrino pulses that were actually 
detected, we have quantitatively addressed the question of which axion masses are consistent 
(or inconsistent) with the experimental data. 

Of all the characteristics of the predicted neutrino pulses, which included total 
number of events, effective neutrino temperature, total energy carried off in neutrinos and 
pulse duration, pulse duration was most sensitively dependent upon the axion mass. In 
particular, for the range of cooling models considered, and an assumed axion mass of 10“ 2 
eV, the predicted pulse duration in the KII detector was less than 1 second and in the IMB 
detector was less than V 2 second — both clearly at variance with the observations. On the 
other hand, for this axion mass, the energy carried off in neutrinos is still -50% or more of 
the total binding energy and the expected number of events were -7-10 for KII and -5 for 
1MB — numbers that are not obviously inconsistent with the actual observations. The reason 
for this is simple: the extended duration of the neutrino burst is connected to the long time 

required for neutrino diffusion to carry the heat trapped in the core of the neutron star to 
the neutrinosphere; with the addition of axion cooling, free-streaming axions from the core 
can rapidly carry off this heat and thereby truncate the late-time part of the neutrino pulse. 

While an axion mass of 10 -2 eV is most certainly ruled out, models that 
incorporate an axion mass of 10" 4 eV are virtually indistinguishable from those without 
axion cooling. From Fig. 4, where At(90%) is plotted vs. axion mass, we see that for all 
cooling models the duration of the neutrino pulses has diminished dramatically for an axion 
mass of about 10~ 3 eV, to less than -6 sec in the KII detector and to less than -2.6 sec in 
the IMB detector. Moreover, for an axion mass of 3x1 0“ 3 eV, the predicted pulse durations 
are less than -2.4 sec (KII) and -1.2 sec (IMB). To summarize then, the neutrino detections 
made by the IMB and KII detectors most emphatically rule out an axion mass of 10” 2 eV, 
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most likely preclude an axion mass as large as 10" 3 eV, and in no way preclude an axion 
mass as small as 10' 4 eV, in agreement with the conclusions reached in Refs. 4 and 5. 21 

We should mention the uncertainties associated with our analysis. First, we have 
relied upon purely theoretical models of the birth and initial cooling of the newly born 
neutron star associated with SN1987A. Since the post-collapse densities are supranuclear, 
there is great uncertainty as to the equation of state. In addition, there is the question of 
the efficiency of the shock at ejecting the outer mantle: how much material eventually 

rains back in on the neutron star? By considering a range of possible post-collapse models, 
we have tried to account for our ignorance, and, as Figs. 2-4 demonstrate, our conclusions 
are very robust and insensitive to the detailed collapse model. For all the models 
considered, the neutrino burst duration decreases dramatically around 10“ 3 eV. 

Perhaps more important are the uncertainties associated with the axion emission 
rate itself. Within the assumptions (the one-pion exchange approximation (OPE)), the 
emission rate has been computed quite accurately (to better than 20%). However, one must 
question the validity of the OPE approximation in general: diagrams involving two-pion and 
other meson exchange may be important. 22 In addition, at supranuclear densities, collective 
nuclear effects, pion condensates, or quark matter in the core might modify the axion 
emission rate and the EOS. It is somewhat reassuring, however, that the axion emission rate 
(<5 a ) is proportional to the axion mass squared. This then means that a factor-of-10 error in 
calculating c? a translates into only a factor-of-3 error in any quoted axion mass limit. 

In sum, in the free-streaming regime (axion masses ^ 0.02 eV), we have shown 
that the existence of an axion more massive than -10" 3 eV would have resulted in neutrino 
pulses of unacceptably short duration in both the KII and IMB detectors. Because the axion- 
nucleon coupling is largely insensitive to whether the axion is hadronic or DFS, this result 
holds for both types of axions. For the hadronic axion, this improves the present mass 

constraint by some three orders of magnitude or more, while for the DFS axion, the 
improvement is only ~1 order of magnitude. Due to axion trapping, an axion more massive 
than ~2 eV may not be precluded by the neutrino burst observations; since other 


astrophysical arguments preclude a DFS axion of this mass, this is only relevant to the 
hadronic axion. Work to address the trapped regime (m a ^ 0.02 eV) is in progress. 
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Table 1. 

Summary of collapse models with 

■ 19 - 

axion emission 

r"' : " J - T 

C: rO’J 

'' 4 \ J--4 

A :,-AtdTY 


m a (eV) 

Number of events 
expected 

KI1 1MB 

Emitted energy 
(10 51 ergs) 
Axions Neutrinos 

At(90%) 

(seconds) 

KII IMB 

Duration of 
calculations 
(seconds) 



Model A: 

(1.3 - 

1.8 Mq, stiff) 




0.0 

15.28 

6.74 

0.0 

328.1 

9.0 

4.0 

20 

10- 4 

15.09 

6.67 

0.1 

324.8 

9.0 

4.0 

20 

3x 1 0 -4 

14.89 

6.63 

14.7 

319.3 

8.0 

3.5 

20 

io- 3 

13.00 

6.10 

77.0 

277.0 

4.0 

1.6 

20 

3xlO" 3 

11.30 

5.51 

121.6 

238.3 

1.6 

1.0 

20 

IO" 2 

9.59 

4.62 

147.2 

215.5 

1.0 

0.5 

14 



Model B: 

(1.3 - 1 

•5 M a stiff) 




0.0 

11.16 

5.45 

0.0 

228.4 

9.5 

4.5 

20 

io- 4 

11.11 

5.37 

2.6 

226.8 

9.5 

4.3 

20 

3xl0 -4 

10.95 

5.37 

7.6 

224.0 

9.2 

4.1 

20 

io- 3 

9.65 

4.97 

45.3 

195.2 

6.0 

2.6 

20 

3xlO -3 

7.73 

4.35 

94.9 

150.5 

2.4 

1.2 

20 

IO" 2 

6.17 

3.96 

127.6 

118.5 

1.0 

0.6 

18 



Model C: 

(1.3 - 1 

■5 Mq, soft) 




0.0 

11.63 

5.84 

0.0 

229.8 

11.0 

5.7 

20 

io- 4 

11.62 

5.84 

7.2 

229.2 

11.0 

5.7 

20 

3xl0- 4 

11.03 

5.64 

36.8 

217.0 

9.1 

4.8 

20 

IO" 3 

9.13 

5.02 

99.6 

176.3 

4.3 

2.0 

20 

3xI0 -3 

7.53 

4.44 

139.6 

141.5 

1.8 

1.0 

20 

IO" 2 

6.40 

3.82 

159.1 

122.2 

1.0 

0.5 

10 
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FIGURE 

Fig. I. 
Fig. 2. 

Fig. 3. 

Fig. 4. 


Fig. 5. 


Fig. 6. 


CAPTIONS 

Feynman diagrams for nucleon-nucleon, axion bremsstrahlung. 

The total neutrino energy (E y ) and the total axion energy (E a ) lost (after 20 
seconds) vs. axion mass, m a , for Models A (solid), B (dashed), and C (dotted). The 
energies are in units of 10 53 ergs and m a is in eV. 

The total expected number of capture events (after 20 seconds) in the IMB 
detector (Nj) and the KIl (N K ) vs. the axion mass, in eV. Models A, B, and C are 
the solid, dashed, and dotted curves, respectively. 

The time required to accumulate 90% of the total number of expected u e capture 
events, At(90%), in seconds, vs. the axion mass, m a , in eV for the IMB and the KII 
detectors for models A (solid), B (dashed), and C (dotted). Note the precipitous 
drop in At(90%) at m a ~ 10 -3 eV. 

Snapshots of the matter temperature (T) in MeV vs. enclosed baryon mass (M) in 
solar masses for Model A without axion emission. The initial model (t - 0) is the 
bottom curve. The snapshots are every 100 milliseconds for the first 2 seconds 
and then every 2 seconds until the end (20 seconds). The compression spike can 
be seen to first grow, then diffuse into the center, and finally begin to decay after 
most of this energy has diffused to the neutrinosphere and is radiated away. 

Same as Fig. 5, but for m a = 10 -2 eV. Note that axion cooling is so effective that 
the inner core never heats up, and thus the energy reservoir which should power 


the late time neutrino emission does not exist. 
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